INVESTIGATION 



Multiple Mechanisms Contribute to Lateral Transfer 
of an Organophosphate Degradation (opd) Island 
in Sphingobium fuliginis ATCC 27551 

Emmanuel Vijay Paul Pandeeti, 1 Toshisangba Longkumer, 1 Deviprasanna Chakka, 
Venkateswar Reddy Muthyala, Sunil Parthasarathy, Anil Kumar Madugundu, Sujana Ghanta, 

Srikanth Reddy Medipally, Surat Chameli Pantula, Harshita Yekkala, and Dayananda Siddavattam 2 

Department of Animal Science, School of Life Sciences, University of Hyderabad, Hyderabad, 500046, India 



ABSTRACT The complete sequence of pPDL2 (37,317 bp), an indigenous plasmid of Sphingobium fuliginis KEYWORDS 

ATCC 27551 that encodes genes for organophosphate degradation (opd), revealed the existence of a site- integrative 
specific integrase (int) gene with an attachment site attP, typically seen in integrative mobilizable elements conjugative 
(IME). In agreement with this sequence information, site-specific recombination was observed between elements (ICE) 

pPDL2 and an artificial plasmid having a temperature-sensitive replicon and a cloned attS site at the 3' genomic islands 

end of the seryl tRNA gene of Sphingobium japonicum. The opd gene cluster on pPDL2 was found to be catabolic 
part of an active catabolic transposon with mobile elements y4qE and Tn3 at its flanking ends. Besides the transposons 

previously reported opd cluster, this transposon contains genes coding for protocatechuate dioxygenase phosphotrieste- 
and for two transport proteins from the major facilitator family that are predicted to be involved in transport rase (PTE) 

and metabolism of aromatic compounds. A pPDL2 derivative, pPDL2-K, was horizontally transferred into organophos- 
Escherichia coli and Acinetobacter strains, suggesting that the oriT identified in pPDL2 is functional. A well- phates 
defined replicative origin (or/V), repA was identified along with a plasmid addiction module re/B/re/E that 
would support stable maintenance of pPDL2 in Sphingobium fuliginis ATCC 27551. However, if pPDL2 is 
laterally transferred into hosts that do not support its replication, the opd cluster appears to integrate into 
the host chromosome, either through transposition or through site-specific integration. The data presented 
in this study help to explain the existence of identical opd genes among soil bacteria. 



Microbial metabolism of organophosphates (OP) attracted the atten- 
tion of microbiologists as it contributes to the elimination of toxic OP 
insecticide residues from agricultural soils (Singh 2009). Microbial 
enzymes involved in degradation of OP compounds are divided into 
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three major groups: organophosphate acid anhydrolases (OAA), 
phosphotriesteases, and methyl parathion hydrolases (Singh 2009). 
The OAAs were later shown to be prolidases involved in hydrolysis 
of the peptide bond of dipeptides with proline at their C terminus 
(Cheng et al 1997; DeFrank and Cheng 1991; Vyas et al. 2010). Nerve 
agents like soman (GD; O-pinacolyl methylphosphonofluoridate), 
sarin (GB; O-isopropyl methylphosphonofluoridate), GF (O-cyclohex- 
ylmethylphosphonofluoridate), and tabun (GA; ethyl JV,N-dimethyl- 
phosphoramidocyanidate) (Cheng et al. 1999) have been shown to be 
fortuitous substrates of prolidases (Cheng et al. 1999). Unlike the 
prolidases, no physiological substrate has yet been found for the other 
two OP hydrolyzing enzymes. Phosphotriesteases (PTE), encoded by 
the organophosphate degradation (opd) gene, are membrane-associ- 
ated metalloenzymes and contain divalent Zn ions at their catalytic 
site (Benning et al. 2001; Mulbry and Karns 1989). The PTEs catalyze 
hydrolysis of the ester linkage found in structurally diverse groups of 
organophosphates, including nerve agents (Benning et al. 1994; Cheng 
and DeFrank 2000; Cho et al. 2004; Lai et al. 1995). The PTE purified 
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Table 1 Bacterial strains and plasmids 



Strain or Plasmid 



Genotype or Phenotype 



Reference or Source 



£ coli DH5a 

E coli BL21 DE3 

£ coli EC100D p/r+1 16 

£ coli HB101 
£ co/i'S17-1 

Sphingobium fuliginis ATCC 27551 

Ac/netobacter sp.DS002 

pRSET-A 

pBluescript KS(II) 

pET-23b(+) 

pMMB206 

pJQ210 

pKD46 

pRK2013 

pSM3 

pPDL2 

pPDL2-K 

pPDL2-KT 
pP4l 
pSFT1 
pSFT2 

pSDP7 
pSDP8 
pSDP9 



supE44zl/acU169 (080 lacZAMIS) hsdRM recA-\ endA-\ 

gyrA96 th/1 re/A1 
hsdS gal(\clts857 indl sam7 ninSlac UV5 17 gene 1 
F mcrA A(mrrhsdRMSmcrBC) (p80dlacZAMT5 AlacX74 recAl 

endA1 araD139 A(ara, /eu)7697 galU galK\- rpsL (StrR) nupG 

pir-1 16(DHFRJ 

P mcrB mrr hsdS20fr E ' m B j recA13 leuB6 ara-14 proA2 lacYt 

galK2 xyl-5 mtl-1 rpsL20(Sm R ) glnV44 A" 
thi pro hsdR hsdM recA RP4 2-Tc::Mu-Kn r ::Tn7(Tp r Sp r Sm r ) 
Wild type strain, Sm r , Pm r , OPH + 
Cm r , Sm r , Ben + 
Amp r , expression vector 
Amp r , /acZ + , cloning vector 
Amp r , expression vector 

Cm r , a low copy number broad host range expression vector 
Gm r , sacB + 

Amp r , Red recombinase expression plasmid 

Kmr, ColE1 -based restricted-host-range helper plasmid 

Amp r , Tc r , opd gene replaced with opd::tet in plasmid pSM2 

opd + , 37.3Kb indigenous plasmid 

opd + - Km r , pPDL2 having single minitransposon 

EZ-Tn5<R6K-yor//KAN-2> insertion 
opd*- Km r , Tet r , pPDL2-K with replacement of opd with opd::tet 
oriT sequence of pPDL2 cloned in pBluescript KS(II) 
repB gene cloned in pRSET-A as BamHI and Xhol fragment 
Tc r , Mini replicon generated by ligating oriV and repA region 

of pPDL2 to opdy.tet cassette 
Amp r , temperature sensitive cloning vector 
Amp r , attB sequence cloned in pSDP7 as H/ndlll fragment 
ligA and ligB genes cloned in pET-23b(+) as Ndel and 

H/ndlll fragment 



Hanahan (1983) 

Studier and Moffatt (1986) 
Epicentre Biotechnologies, USA 



Boyer and Roulland-Dussoix (1969) 

Simon et a/. (1983) 
Kawahara eta/. (2010) 
Unpublished results from our lab 
Invitrogen, USA 
Fermentas, USA 
Novagen, USA 
Morales eta/. (1991) 
Quandt and Hynes (1993) 
Datsenko and Wanner (2000) 
Ditta et a/. (1 980) 
Siddavattam et a/. (2003) 
Mulbry eta/. (1986) 
This study 

This study 
This study 
This study 
This study 

This study 
This study 
This study 



from Brevundimonas diminuta hydrolyses parathion, the model OP 
compound, at a rate close to its diffusion limits and is considered to be 
an end point in enzyme evolution (Caldwell and Raushel 1991; Scan- 
Ian and Reid 1995). The OPs were recently (65 years ago) introduced 
into agriculture pest management, mainly as replacements for more 
persistent organochlorides (Davies et al. 1985). Therefore, evolution of 
PTEs in such a short period has become an interesting model in which 
to study molecular evolution (Elias and Tawfik 2012; Scanlan and 
Reid 1995). The studies conducted to date have shown quorum- 
quenching lactonases as progenitors of PTEs (Afriat et al. 2006; Elias 
and Tawfik 2012; Roodveldt and Tawfik 2005). 

The third prominent organophosphate-degrading enzymes are 
methyl parathion hydrolases (MPH), encoded by the methyl para- 
thion degradation (mpd) gene. The MPH enzymes known to date 
have been purified from Pseudomonas strains isolated from OP-pol- 
luted Chinese agricultural soils (Liu et al. 2005; Zhang et al. 2006; 
Zhongli et al. 2001). Despite having functional similarity, MPH 
enzymes share no homology with PTEs, indicating the existence of 
structurally independent organophosphate-degrading enzymes among 
soil bacteria (Dong et al. 2005). Unlike organophosphate hydrolase 
(OPH), the MPH has been shown to be a structural homolog of 
^-lactamases, the enzymes that confer resistance to f5-lactam-derived 
antibiotics (Tian et al. 2008). Interestingly these two structurally dif- 
ferent enzymes have an identical active site structure (Dong et al. 
2005). This type of functional convergence pointed to the existence 
of independent paths in the evolution of organophosphate-hydrolyz- 
ing enzymes. Lateral transfer of mpd genes became evident with the 
discovery of identical mpd genes among bacterial strains isolated from 



OP-polluted Chinese soil samples (Zhang et al. 2006), and the mpd 
genes were shown to be part of an active transposon (Wei et al. 2009). 

Identical opd genes have been found in soil bacteria isolated from 
diverse geographical regions. The dissimilar opd plasmids pCMSl and 
pPDL2 were isolated from B. diminuta (pCMSl) and Flavobacterium 
sp. ATCC 27551 (pPDL2) that were enriched, respectively, from soil 
samples collected from agricultural fields of Texas, USA, and the In- 
ternational Rice Research Institute (IRRI), Philippines. They were 
shown to contain identical opd genes (Mulbry et al. 1986; Pandeeti 
et al. 2011), and a 7 kb region around the opd gene apparently con- 
stituted the only identity between these two dissimilar plasmids 
(Mulbry et al. 1986; Pandeeti et al. 2011; Siddavattam et al. 2003). 
However the sequence of the opd gene cluster in the self-transmissible 
pCMSl showed no features of a transposable element (Pandeeti et al. 
2011), and although a transposon-like organization was found in 
pPDL2, the opd cluster of pPDL2 has not been shown to be an active 
transposon (Siddavattam et al. 2003). 

In the present study, we report the complete sequence of pPDL2, 
and we show experimentally that pPDL2 is a mobilizable plasmid 
having the capability to spread the opd cluster laterally by functioning 
both as an integrative mobilizable element and as an active transposon. 

MATERIALS AND METHODS 
Strains and plasmids 

Bacterial strains and plasmids used in the present study are shown in 
Table 1. Oligonucleotide primers used while performing PCR reac- 
tions are listed in supporting information, Table SI. E. coli and Aci- 
netobacter sp. DS002 cells were grown in LB medium at 37° and 30°, 
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respectively. The cultures of Sphingobium fuliginis ATCC 27551 were 
grown in modified Wakimoto medium at 30°. When necessary, the 
antibiotics ampicillin (100 u,g/ml), tetracyclin (30 jjug/ml), gentamycin 
(20 |xg/ml), and kanamycin (30 u-g/ml) supplemented the culture 
medium. The LB sucrose plates were prepared by adding 5% sucrose 
to LB agar plates. Acinetobacter sp. DS002 strains were grown in M9 
medium containing filter-sterilized benzoate (5 mM) as a carbon 
source. Exconjugants of Acinetobacter sp. DS002 strains were selected 
on M9 agar plates supplemented with kanamycin. 

Isolation and rescue cloning of pPDL2 from 
Sphingobium fuliginis ATCC 27551 

Large indigenous plasmids of Sphingobium fuliginis ATCC 27551 were 
isolated by following the Currier Nester protocol with minor modifica- 
tions (Currier and Nester 1976; Pandeeti et al. 2011). The plasmid 
preparation made from Sphingobium fuliginis ATCC 27551 was directly 
used for tagging with minitransposon EZ-Tn5<R6K7ori/KAN-2> us- 
ing the EZTn5 <R6K7ori/KAN-2> insertion kit (Epicenter Biotech- 
nologies, USA) following the manufacturer's protocols. The isolated 
plasmid preparation and minitransposon was taken in equimolar con- 
centrations and incubated with transposase for 2 hr at 37° to complete 
in vitro transposition. After 2 hr, transposase was inactivated by adding 
1 |jl1 of 1% SDS before incubating the reaction mixture for 10 min at 
70°. The reaction mixture was then electroporated into E coli EC100D 
pir-116 cells by setting the electroporator (Gene Pulser, Bio-Rad Labo- 
ratories, USA) at 2.5 kV, 200fl for 4.5 msec. After electroporation, the 
cells were plated on LB plates supplemented with kanamycin to select 
transformants. Colonies having plasmid pPDL2 were identified by per- 
forming colony PCR using opd specific primers. The colonies that gave 
amplification of the opd gene were used to isolate plasmids having 
the minitransposon Tn5<R6K7ori/KAN-2>. The derivatives of plas- 
mid pPDL2, named pPDL2-K, were then isolated from E. coli pir-116 
cells using the BAC isolation protocol (Sambrook et al. 1989). The 
restriction profile of pPDL2-K was generated by digesting with Smal, 
BamHl, Pstl, Sail, Xhol, Hindlll, and EcoRI, and it was compared with 
the restriction profile of pPDL2 (Mulbry et al. 1987) to identify the exact 
point of minitransposon insertion. Further, all fragments of pPDL2-K 
obtained after digestion with Pstl and EcoRI were ligated independently 
into pBluescript (II) KS (+) digested with similar enzymes. 

Generation of the complete pPDL2 sequence 

DNA sequencing was performed by using BigDye Terminator vl.l 
Cycle Sequencing Kit (Applied Biosystems, USA), and an ABI PRISM 
3100 Genetic Analyzer was used following manufacturer's protocol. 
Either pPDL2 derivatives or subclones generated by cloning EcoRI 
and Pstl fragments were used as templates. Sequence reactions were 
generated using transposon-specific sequence primers when pPDL2 
derivatives were used as templates, whereas universal forward and re- 
verse primers were used to generate sequence from subclones. To fill 
gaps, a primer- walking strategy was employed. If the cloned fragment 
size was large, it was digested with different restriction enzymes and 
subcloned into pBluescript (II) KS(+) and were used for sequencing. 

Sequence assembly and annotation of pPDL2 

All sequences were viewed and edited to remove vector sequences by 
using Chromas 2.13 software (www.technelysium.com.au/chromas). 
Sequences were assembled into contigs by using the program Contig 
Express of VectorNTI software (Invitrogen Technologies, USA). The 
assembled sequence of plasmid pPDL2 was annotated using Artemis 
sequence annotation tool (http://www.sanger.ac.uk/resources/software/ 



artemis/) (Rutherford et al. 2000). The open reading frames (ORFs) 
were identified by using the built-in tool of the Artemis software, and 
the start codons in the predicted ORFs were fixed with the help of 
BLAST searches. BLAST searches were made against the nonredun- 
dant database of NCBI using BLAST program (www.ncbi.nlm.nih.gov/ 
BLAST). IS elements, transposons, and their repeat elements were 
identified by doing a pairwise alignment using BlastN program against 
the ISfinder database (http://www-is.biotoul.fr/). 

Prediction of origin of replication (oriV) 

The oriV was predicted based on sequence homology to other anno- 
tated or predicted oriVs and by performing GC skew analysis (Gri- 
goriev 1998). While calculating the GC skew, the ratio of (G-C) to (G 
+C) was calculated per each window of equal length split over the total 
sequence of pPDL2. The GC skew was multiplied by w/c to find its 
dependence on subsequence length, where w and c are lengths of the 
subsequence and total sequence of pPDL2, respectively. As maximum 
and minimum value of GC skew is usually associated with termination 
and origin of replication (Frank and Lobry 1999), a similar approach 
was followed to predict oriV of plasmid pPDL2 origin. 

Prediction of attP sites 

For predicting attP sites, two independent approaches were used. In 
the first approach, sequences of all available plasmids having identical 
integrases were collected from the NCBI database. From bacteria 
having these plasmids, a dataset of all tRNA sequences along with 
their upstream and downstream sequences was created. Pairwise 
alignments were made between the plasmids and the tRNA sequences 
to identify the att sites. 

In a second approach, the genome sequences of bacteria having 
plasmids with integrase genes were collected from NCBI database. 
These genomic sequences were then used to predict genomic islands 
(GI) using Islandviewer software (http://www.pathogenomics.sfu.ca/ 
islandviewer/query.php). The predicted GI sequences were then used 
to make pairwise alignments with sequences of tRNA genes of Sphin- 
gomonas wittichii and Sphingobium japonicum. Short sequence repeats 
that exactly matched the 3' end of tRNA sequences were taken as 
putative att sites. The predicted att sites with a low E-value were 
considered as potential att? and attB sites. The att sites predicted in 
this manner were then used to align with plasmid pPDL2 sequence 
using BlastN to identify attP homologs. Alternatively, tRNA gene 
sequences taken from S.wittichii and S. japonicum were directly used 
to align with the pPDL2 sequence using BlastN to find short sequen- 
ces that perfectly matched the 3' end of tRNA genes. 

Bacterial conjugation experiments 

Biparental mating experiments were done using E. coli pir-116 
(pPDL2-K) as donor and Acinetobacter sp. DS002 as recipient. The 
triparental mating experiments were performed by including E. coli 
strain HB101 (pRK2013) as a helper. Biparental and triparental experi- 
ments were performed following standard protocols described else- 
where (Figurski and Helinski 1979). The exconjugants were selected 
by plating on selection plates having minimal media supplemented 
with 5 mM benzoate and kanamycin. Parent strains treated in a similar 
manner were plated on selection plates and served as negative controls. 

Cloning of synthetic attB site in a temperature- 
sensitive vector 

A temperature-sensitive plasmid vector was constructed by using 
vectors pMMB206 (Morales et al. 1991) and pKD46 (Datsenko and 
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Wanner 2000) as source plasmids. The temperature-sensitive replica- 
tive origin along with the bla gene were amplified as an EcoRVISmal 
fragment by using primers DSF005 and DSF006 as forward and re- 
verse primers. The multiple cloning site and lac reporter system were 
taken from plasmid pMMB206 by digesting it with EcoRV/Dral. 
These two fragments were then ligated to give plasmid pSDP7. The 
seryl tRNA with an attB sequence was synthesized before cloning it in 
pSDP7. Equal concentrations of the overlapping primers DSF001 and 
DSF002 were placed in a PCR tube in boiling water and left until the 
water reached room temperature. The contents were then brought 
down by briefly spinning, and then 2.5 u,mol dNTP, a unit of pfu, 
and 10 pi lOx pfu buffer were added. The contents were made up to 
100 u,l with MilliQ water before incubating the tube at 72° for 10 min. 
An aliquot of 10 ul of the reaction mixture was then taken as a tem- 
plate to amplify attB by performing PCR using the primer set 
DSF003/DSF004 with Hindlll recognition sites at the 5' end. The 
amplicon was then digested with Hindlll and cloned into pSDP7 to 
give pSDP8. 

In vivo integration assay 

Plasmid pSDP8 was transformed into E. coli pir-116 (pPDL2-K), and 
the resultant cells were grown at 30° for 12 hr in LB medium with 
kanamycin and ampicillin. After incubation, cells were collected from 
1 ml of culture and thoroughly washed before reinoculating them in 
fresh LB medium containing only kanamycin. The cultures were then 
grown for a further period of 12 hr at 37° to inhibit the replication of 
the temperature-sensitive plasmid pSDP8. After incubation, the cul- 
ture was serially diluted, and appropriate dilutions were plated on LB 
agar with kanamycin and ampicillin. Plasmids were isolated from the 
resultant colonies, and the formation of a cointegrate due to site- 
specific recombination between attP and attB was determined by 
performing restriction analysis. The attB/P and attPIB sites were 
identified by DNA sequencing using universal forward and reverse 
primers. 

Transposition assay 

For the transposition assay, we created a pPDL2-K derivative by 
replacing the opd gene with opdv.tet through homologous recombina- 
tion. E. coli pir-116 (pPDL2-K) was first transformed with the tem- 
perature-sensitive plasmid pKD46, which codes for red recombinase, 
in the presence of arabinose (Datsenko and Wanner 2000). The mid 
log-phase cultures obtained by growing in LB medium with 1 mM 
arabinose were made electrocompetent and immediately used to elec- 
troporate with insertionally inactivated opd (opdv.tet) taken as a Pstl 
fragment from pMS3 (Siddavattam et at 2003). The tetracycline-re- 
sistant colonies were then used to perform colony PCR using opd- 
specific primers. The colonies identified as carrying opdv.tet were 
grown on them at 37° to eliminate temperature-sensitive pKD46. 
The ampicilin-sensitive, tetracycline-resistant colonies were then used 
to isolate the pPDL2-K derivative pPDL2-KT and to confirm the 
presence of opdv.tet by restriction analysis. 

The transposition assay used in the present study is described 
elsewhere (Siddavattam et al. 2003), except that the sacB containing 
pJQ210 (Quandt and Hynes 1993) was used as a reporter plasmid and 
plasmid pPDL2-KT served as donor of opd cluster. E. coli (pJQ210) 
cells were made electrocompetent and plasmid pPDL2-KT was elec- 
troporated following standard protocols (Dower et al. 1988). The 
electroporated cells were then grown in LB medium having no anti- 
biotics for 2 hr. The cells were then plated on LB plates having 
gentamycin, sucrose, and tetracyclin. As plasmid pPDL2-KT requires 



a pir background, the tetracycline-resistant colonies found on sucrose 
plates were considered to have been generated due a transposition 
event disrupting the sacB gene. These colonies were then used to 
isolate plasmid and to analyze junction sites by performing restriction 
profile and PCR amplification. The PCR amplification was done by 
using sacB-specific primers (DSF0015/DSF0016) and primers 
(DSF0017/DSF0018) designed using pPDL2 sequence found upstream 
of the predicted Tn3-specific terminal repeats. 

Interactions between or/T and Mob protein 

The predicted repB gene downstream of oriT was amplified using 
a primer set DSF009/DSF0010 appended with BamHl and Xhol sites. 
The amplicon was then cloned in expression vector pRSETA as 
a BamHl/Xhol fragment to generate recombinant plasmid pSFTl. E. 
coli BL21 (pSFTl) cells were induced to express RepB c _ 6His and affin- 
ity purified following procedures described elsewhere (Pandey et al. 
2009). A 250 bp DNA fragment containing the oriT region was am- 
plified using a primer pair (DSF0019 and DSF0020) designed by 
taking the sequence flanking oriT of pPDL2. The PCR amplicon 
containing oriT was end-labeled with [7- 32 P] ATP using T4 polynu- 
cleotide kinase following established protocols (Sambrook et al. 1989). 
While performing mobility shift assay, the pure labeled oriT (2 pmol) 
was taken in 20 pi of binding buffer [20 mM Tris-HCl (pH 8.0), 1.0 
mM EDTA, 6 mM MgCl 2 , 50 mM KC1, 50 fjLg/ml BSA, 5% (w/v) 
glycerol, 5 p.g/ml of herring sperm DNA] and incubated for 20 min at 
25° with various concentrations (0 ng, 100 ng, 250 ng, 500 ng and 
1000 ng) of RepB C -6His- The reaction mixture was then resolved on 
a 5% (w/v) native polyacrylamide gel, and the mobility of labeled oriT 
was captured on X-ray film by performing autoradiography. Control 
reactions were performed either by omitting RepB c _ 6His or by includ- 
ing increased concentrations (0.2 to 2 pmol) of cold oriT. 

Expression of the lig operon 

The lig operon comprising the ligB and HgA genes was amplified from 
plasmid pPDL2 using primer set DSF0013/DSF0014 appended with 
Ndel-Hindlll, respectively. The Ndel site in the forward primer cor- 
responds to the 5' region of ligB and overlaps with its start codon. 
Similarly, the reverse primer designed by taking the 3' end of HgA 
contains a Hindlll site created by changing the stop codon. The 
amplicon containing the ligB A genes was then cloned in pET23b as 
an Ndel/Hindlll fragment to generate plasmid pSDP9. The cloning 
strategy followed to generate pSDP9 places both ligB and HgA under 
the transcriptional and translational signals of the expression vector. 
As the stop codon of HgA was eliminated to facilitate fusion of HgA 
with the vector-coded histidine tag, the LigA alone was expressed with 
C-terminal histidine tag. As active protocatechuate dioxygenease con- 
tained both LigB and LigA subunits, affinity purification using Ni- 
affinity column was expected to yield LigBA complex. Expression of 
LigB A was monitored both by measuring dioxygenease activity (Ono 
et al. 1970) and by detecting LigA c _6His from Western blots using anti- 
his antibody (Pandey et al. 2009). 

Enzyme assays 

Protocatechuate-4,5-dioxygenase assay was performed using purified 
LigBA complex. The reaction mixture contained protocatechate (20 
pJVI) and 5 p,g of purified LigBA complex in 100 mM Tris acetate 
buffer (pH 7.2). Formation of the ring clevage product, 2-hydroxy-4- 
carboxymuconic semialdehyde, was measured spectroscopically fol- 
lowing established procedures (Ono et al. 1970). Similarly, formation 
of 2-hydroxymuconic semialdehyde was monitored at 375 nm when 
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catechol was used as a substrate (Nozaki et al. 1970). The organo- 
phosphate hydrolase (OPH) assay was performed using methyl para- 
thion as substrate. The OPH activity was measured by estimating the 
concentration of p-nitrophenol (Pakala et al. 2007). Specific activity of 
enzymes was expressed as micromoles of product formed per minute 
per milligram of protein. 

RESULTS 

Flavobacterium sp. ATCC 27551 was one of the first organophos- 
phate-degrading bacteria isolated from the agricultural soils of the 
International Rice Research Institute, Manila, Philippines (Sethuna- 
than and Yoshida 1973). Recentiy it was reclassified as Sphingobium 
fuliginis ATCC 27551 (Kawahara et al. 2010). The OP degradation 
property has been shown to be encoded on one of its four indigenous 
plasmids, pPDL2 (Mulbry et al. 1986), the other three not being in- 
volved in this activity. Lateral gene transfer of opd is evident due to the 
existence of identical opd genes in bacterial strains belonging to dif- 
ferent taxonomic groups. However, no systematic studies are available 
to explain the mechanism of lateral transfer of the opd gene cluster. 
Because pPDL2 contains a transposon-like opd cluster, we have taken it 
as a model system to study horizontal gene transfer (HGT) of the opd 
cluster. Initially pPDL2 was isolated from the other indigenous plasmids 
in S. fuliginis ATCC 27551 by performing rescue cloning. After in vitro 
transposition using the minitransposon Tn5 <R6K7ori-Kan2>, 
a pPDL2 derivative (pPDL2-K) having a single minitransposon inser- 
tion was identified by colony PCR using opd gene-specific primers. 

The complete nucleotide sequence of 37,317 bp plasmid pPDL2 is 
deposited in GenBank (accession no. JX31-2671). Its GC content 
(62.37%) is very similar to the GC content of indigenous plasmids 
pCHQl (63%), pUTl (63.7%), and pUT2 (61.0%) found in S. japo- 
nicum (Nagata et al. 2010, 2011). The GC content, GC skew, and the 
length of the plasmid indicates multiple changes in the sequence of 
pPDL2 either due to acquisition of foreign DNA or rearrangement of 
its own sequences (Figure 1). A total of 48 protein-coding regions 
were annotated, which represents a coding density of 91.3% (Table 
S2). Out of the 48 orfs, 18 of them were annotated as ORFs coding for 
hypothetical proteins, and 7 of them code for either transposases or 
associated resolvases (Table S2). Based on comparison with database 
entries, the ORFs coded by plasmid pPDL2 can be divided into a num- 
ber of functional modules. Prominent among them are modules for 
replication and partition, mobilization, integration, and OP degrada- 
tion, as well as mobile genetic elements. 

Replication initiation proteins 

In pPDL2, two ORFs, orf22 and orf45, show high homology to rep- 
lication initiator proteins. The protein coded by orf22 (367aa) shows 
absolute identity to the replication initiation protein RepA of plasmid 
pUTl of S. japonicum pUT26S (Nagata et al. 2010). It contains a se- 
quence that shows similarity with the well-defined RPA and Rep-3 
domains. The proteins having these conserved domains have been 
shown to play a key role in replication initiation of plasmids (Wickner 
et al. 1991). As in plasmid pUTl, in the upstream of repA, an oriV 
sequence was identified with characteristic iterons, an AT-rich se- 
quence, and a dnah box (Figure 2). Iterons play a key role both in 
initiation and regulation of replication in theta- replicating plasmids 
(del Solar et al. 1997; Nordstrom 1990). There is a significant simi- 
larity among the oriV sequences found in plasmids pPDL2, pUTl, and 
pPSlO (Figure 2). If the structural conservation of both oriV and RepA 
is taken into consideration, they qualify to have all the potential to 
replicate plasmid pPDL2 (Figure 2). To gain further evidence on their 
role in pPDL2 replication, we generated a mini-replicon by ligating 



the rep A oriV region (18,815-20,914) to a tetracycline-resistant gene. 
When electroporated into Sphingobium fuliginis ATCC 27551, the 
resultant mini-replicon (pSFT2) gave tetracycline-resistant colonies, 
indicating the predicted oriV and repA contribution for the replication 
of pPDL2. Further, the generated mini-replicon remained in the cells 
even after growing them in the absence of antibiotic for several gen- 
erations. Existence of the mini-replicon was also evident in plasmid 
preparations made from the transformed cells (data not shown). The 
RepA sequences of theta-replicating plasmids have been showed to 
follow a characteristic phylogenetic pattern and, to some extent, serve 
as signatures when assigning incompatibility groups to the plasmids 
(del Solar et al. 1997; Suzuki et al. 2010). Therefore, the RepA se- 
quence of plasmid pPDL2 was aligned with the available RepA 
sequences in the GenBank (Suzuki et al. 2010). RepA of pPDL2 has 
more than 50% identity covering the complete sequence of RepA 
sequences of the IncP group of plasmids, pADP-1, pEST4011, and 
pUBl (Suzuki et al. 2010). Such similarity is not seen with RepA 
sequences of other incompatibility groups of plasmids, suggesting that 
pPDL2 could be an IncP group of plasmid. 

The protein coded by orf24 of pPDL2 has 100% sequence coverage 
and identity with the partitioning protein Par A of plasmid pUTl of S. 
japonicum UT26S. ParA, in association with ParB, plays a major role 
in plasmid segregation. However, no parB homolog is found in 
pPDL2. In the absence of parB, pPDL2 may depend on a ParB protein 
encoded either by the chromosome or by other indigenous plasmids 
present in S. fuliginis ATCC 27551. Interestingly, a well-defined plas- 
mid addiction module was identified in pPDL2. The ORFs orf29 and 
orf30 code for proteins showing more than 95% identity to the RelB/ 
RelE toxin/antitoxin system found in a variety of Gram-negative bac- 
teria and archaea (Gotfredsen and Gerdes 1998; Gronlund and Gerdes 
1999). The RelE and RelB pair of E. coli is well characterized. The 
1 lkD RelE protein acts as toxin to bacterial cells, and a similarly sized 
relB product counteracts the toxic effect of RelE (Gotfredsen and 
Gerdes 1998; Gronlund and Gerdes 1999). Presence of partition pro- 
tein ParA, along with a post-segregational killing mechanism (RelE/ 
RelB), suggests faithful partitioning and maintenance of pPLD2 in 
daughter cells of Sphingobium fuliginis ATCC 27551. 

Plasmid pPDL2 is a mobilizable plasmid 

Since oriT is an essential element for mediating horizontal gene trans- 
fer we searched for a similar sequence in pPDL2. A number of func- 
tionally validated oriT sequences have been shown to contain typical 
secondary structures in a characteristic GC rich region (Lee and 
Grossman 2007) and such an onT-like sequence (34823-35011) was 
identified in pPDL2. The predicted oriT has very high similarity to an 
oriT sequence of a self-transmissible plasmid pLBl of S. japonicum 
UT26S (Miyazaki et al. 2006). If the location of oriT is taken into 
consideration, along with the presence of four tandem repeats imme- 
diately downstream, the predicted sequence has all features of a func- 
tional oriT (Figure 3A). Interestingly the putative oriT sequence is 
found upstream of an ORF (orf45) that codes for a protein showing 
high similarity (98%) to the plasmid-borne RepB sequences of Sphin- 
gobiaceae members (Bertalan et al. 2009). 

The predicted RepB sequence includes a Rep-3 domain, involved 
in typical nicking-closing like activity (Dyda and Hickman 2003; 
Moscoso et al. 1995). Proteins having a Rep-3 domain also perform 
strand-transfer reactions generally seen during rolling circle (RC) rep- 
lication (Moscoso et al. 1995; Ruiz-Maso et al. 2007). RC replication is 
also seen during horizontal transfer of indigenous plasmids. However, 
it is initiated from a well-defined oriT sequence (Pansegrau and Lanka 
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Figure 1 Physical map of plasmid pPDL2. Outer and inner circles indicate proteins encoded by sense and anti-sense strands, respectively. Third 
circle indicates mobile elements and repeat sequences. Direct (DR) and Inverted (IR) repeats are shown with filled red and green triangles, 
respectively. Tn3-specific repeats appear with filled purple triangles in the fourth circle. The fifth circle shows GC content across the plasmid 
sequence. The sixth and seventh circles represent GC-skew in sense and anti-sense strands. 



1991). Since we have identified both oriT and RepB in the sequence of 
pPDL2 we have attempted to validate their involvement in the lateral 
transfer of pPDL2. We initially assessed whether pPDL2 is a mobiliz- 
able plasmid by performing triparental mating experiments, involving 
E. coli pir-116 (pPDL2-K) as donor, with Acinetobacter sp. DS002 
and E. coli (pRK2013) as recipient and helper strains respectively. 
Kanamycin-resistent exconjugants of Acinetobacter sp. DS002 have 
appeared on a selection plate at a frequency of 3.4xl0~ 5 per recipient. 
As expected, all exconjugants carried the opd gene and showed OPH 
activity, properties not found in Acinetobacter sp. DS002 (Figure 3D). 



Mobility of plasmid pPDL2 was also observed when E. coli BL21 was 
used as recipient. Though pPDL2 is stably replicated in Acinetobacter 
sp. DS002, it failed to replicate in an E. coli pir negative background. 
Although kanamycin resistant E. coli colonies were found, they sub- 
sequently lost the plasmid. In certain colonies that retained the ability 
to grow on kanamycin we could not detect plasmid from these cells 
suggesting that either complete pPDL2-K or the region containing 
kanamycin gene has probably integrated into the genome. 

To understand RepB and oriT interactions better we performed 
mobility shift assays by incubating the PCR amplicon carrying oriT 
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with purified RepB C -6His (Figure 3B). There was a RepB c _ 6His depen- 
dent shift in the electrophoretic mobility of oriT (Figure 3C) and no 
shift in presence of a BSA control, indicating specific interactions 
between RepB C - 6 ms and oriT. Supporting this observation the inten- 
sity of shifted oriT band got proportionately reduced when increased 
concentrations of cold oriT was added to the reaction mixture (data 
not shown). To assess if RepB might function as a relaxase (Mob) the 
sequence of RepB was aligned with other relaxase sequences (Acht- 
man 2011; Francia et al 2004; Garcillan-Barcia et al. 2009; Gugliel- 
mini et al. 201 1) and a weak homology was observed with relaxase of 
the MOB F plasmid. The role of RepB in lateral transfer of pPDL2, 
therefore needs to be further investigated. Despite not detecting tra 
genes in the sequence of pPDL2, experiments were also performed to 
assess whether pPDL2 is self-transmissible but no exconjugants were 
observed from a biparental mating using E. colipir-116 (pPDL2-K) as 
donor and Adnetobacter sp. DS002 as recipient. Hence, in summary, 
our data clearly indicate that pPDL2 is a mobilizable but not a self- 
transmissible plasmid. 

Integration module 

An interesting observation in the sequence of pPDL2 is the existence 
of an uniquely organized integrase module (nucleotide position 26697 
to 30249). The module contains two units, CIP-I and CIP-II, each 
consisting of copG, int and pgm genes and having divergent transcrip- 
tion orientation (Figure 4A). The integrase of CIP-I is a 328 amino 
acid protein with strong identity (86%) to a plasmid pUTl -borne 
phage integrase of S. japonicum UT26 and 70% identity to a site- 
specific recombinase of Pseudomonas syringae DC 3000. The CopG/ 
MetJ/Arc family regulatory protein (subsequently referred to as 
CopG) coded by pPDL2 shows highest similarity to the CopG protein 
coded by plasmid pUTl of S. japonicum UT26S. The 131 residue 
protein shows homology to a number of other proteins belonging 



to CopG/MetJ/Arc family of transcriptional regulators that can act 
both as transcriptional repressors and activators (del Solar et al. 
1995; del Solar and Espinosa 1992; Gomis-Ruth et al 1998; He 
et al. 2002). In a well-characterized Streptococcal plasmid, pMV185, 
CopG has been shown to act as a repressor of repB by binding to 
a pseudosymmetirc region overlapping the -35 sequence of the cr 70 - 
dependent promoter, thereby preventing transcription of the copG, 
repB operon in plasmid pMV158 (Costa et al. 2001). 

In pPDL2, the copG and int genes appear to be organized as one 
operon in both CIP-I and CIP-II units. A potential a 70 -dependent 
promoter was only seen upstream of each of the copG sequences, and 
the stop codon of copG overlaps with the start codon of the int genes 
suggesting that these two genes are cotranscriptional. Due to the 
existence of opposite transcriptional orientation, an overlap of 5bp 
is seen between the putative copGl and copG2 promoters (Figure 
4A). As CopGl proteins are known to act as transcription repressors, 
the intergenic region found between the copGl, intl and copG2, intl 
operons was analyzed to identify cis-elements that can serve as target 
sites for CopG proteins (Figure 4A). In studies conducted by del Solar 
and his associates, CopG is shown to bind to a pseudosymmetric cis- 
element present overlapping -35 region (del Solar et al. 2002). Inverted 
repeats are also shown to be putative binding sites of CopG in Sulfo- 
lobus neozealandicus (Greve et al. 2005). Interestingly, in both CIP-I 
and CIP-II units, typical CopG binding sites were predicted upstream 
of the start codon of copG and putative promoter sequence (Figure 
4A). Transcriptional organization of copG and int genes and the ex- 
istence of such CopG binding sites, clearly indicates involvement of 
CopG in regulation of expression of the integrase in S. fuliginisATCC 
27551. 

In addition to these two transcriptionally coupled integrase and 
copG genes an additional ORF that codes for a protein showing iden- 
tity (98%) to phosphoglycerate mutase of S. japonicum UT26S is 
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present at the 3' end of the int gene in each of the two units of the 
integrase modules (Figure 4A). Though, the phosphoglycerate mutase 
gene in CIP-I module is found 1.1 Kb away from the stop codon of 
integrase, in CIP-II, it is found immediately downstream of transcrip- 
tion terminator sequence of intl gene. This 206 residue protein shows 
high identity (98%) to phosphoglycerate mutase encoded on plasmid 
pUTl of S. japonicum UT26S. Phosphoglycerate mutases (PGMs) are 
very well characterized group of enzymes found both in prokaryotes 
and eukaryotes. PGMs are basically transferases involved in transfer of 
phosphate group (Parkinson and Kofoid 1992). Conversion of glyc- 
eraldehyde 3-phosphate to glyceraldehyde 2-phosphate, by transfer of 
phosphates from the 3 rd position to the 2 nd position, is the classical 
biochemical reaction catalyzed by this group of enzymes. However, 
PGMs containing an HP-PGM-like domain are also known to be 
involved in signal transduction process (Matsubara and Mizuno 
2000). The histidine present in the catalytic site of this group of 
enzymes undergoes phosphorylation during the signal relay process 
(Parkinson and Kofoid 1992). The strong link between PGM, CopG 
and integrase suggests a related function of these proteins, and a sim- 
ilar organization is also seen in S. japonicum UT26S plasmids pCHQl 
(NCBI Reference no: NC014007) and pUTl (NCBI Reference no: 



NC014005), and Sphingomonas wittichi RW1 (pSWITOl) adding 
strength to the proposed hypothesis. 

To assess whether the CIP-I and CIP-II units were generated by gene 
duplication we compared the sequences of these two units. If they are 
products of gene duplication there should have been absolute identity 
between the proteins coded by these two units of the integrase module. 
However there is only 67% identity between these two proteins. The 
N-terminal region of CopG2 encoded by CIP-II unit was found to be 
much longer. In addition to the size difference, there was no sequence 
conservation at the C-terminus of the protein. Nevertheless, the central 
regions of CopGl and CopG2 are almost identical. Such diversity in the 
primary sequence suggests a divergent origin of CopGl and CopG2. 
A similar situation was seen when intl and pgml were aligned with intl 
and pgm2. Between intl and intl only 80% identity was observed al- 
though this identity continued throughout the sequence. The percent 
identity between PGM1 and PGM2 is high (91%), the difference only 
being seen at the C-terminus and N-terminus. Considering the diversity 
in the primary sequence of the proteins coded by CIP-I and CIP-II of 
the integrase module, the plasmid pPDL2 is proposed to have acquired 
the CIP-I and CIP-II integrase modules from independent sources, 
possibly through a unique recombination process. 
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Plasmid pPDL2 is an integrative mobilizable element 

Having identified an integrase module, we performed experiments to 
determine whether site-specific integration of pPDL2 takes place at an 
artificially created attB site cloned in temperature- sensitive plasmid 
pDSP8. Initially we performed bioinformatic studies to identify attP 
sites on the plasmid pPDL2. While perfoming these predictions, we 
took tRNA genes from the genome sequence of S. japonicum UT26S, 
whose total genome sequence is known (Nagata et at 2010). The 
reasons for taking tRNA genes from S. japonicum are that parts of 
pPDL2 show considerable homology with its indigenous plasmid 
pUTl and that it encodes an integrase module very similar to the 
one present in pPDL2. These tRNA gene sequences were then blasted 
with the total sequence of pPDL2 to obtain a library of identical 
sequences with a minimum length of 10 bp. Among these sequences, 
the ones found in intergenic regions were predicted to be potential 
attP sites. The seryl tRNA gene of S. japonicum has sequence identity 
with the predicted attP at its 3' end and served as an attB sequence 
[Figure 4B, part (i)]. E. coli pir-116 strains carrying pPDL2-K and 
pSDP8 were grown independentiy for 12 hr at 30° and 37°, and then 



plated on kanamycin and ampicillin plates. A lawn of cells were seen 
on plates seeded with the 30°-grown culture, whereas only a few 
colonies appeared in plates seeded with the 37°-grown culture. In 
the absence of attB, no colonies appeared, suggesting that the gener- 
ated colonies were due to occurrence of site-specific recombination 
between attB and attP sites [Figure 4B, part (ii)]. The restriction pro- 
file of plamids isolated from randomly picked colonies showed in- 
tegration of pSDP8 in a 12.0 kb Pstl fragment of pPDL2 (Figure 4, C 
and D), and sequence analysis revealed attPBL and attBPR sequences, 
suggesting that the attP site predicted at the 3' end of intl served as an 
attP site in pPDL2. 

Degradation module 

In our previous studies, we analyzed the 7 kb DNA region of pPDL2 
that shows sequence identity with the pCMSl -borne opd cluster from 
B. diminuta, and we identified a transposon-like organization of the 
opd gene cluster in S.fuliginis ATCC 27551 (Siddavattam et al. 2003). 
Such structural organization and the existence of identical opd clusters 
in bacterial strains representing diverse taxonomic groups suggested 
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horizontal transfer of the opd cluster through transposition (Figure 
SI). However, in our previous studies, we were unable to detect any 
transposition events, and we concluded that it is a defective trans- 
poson due to the existence of a truncated tnpA gene in the Tn3 
element identified at the flanking end of the opd cluster. 

The full sequence of pPDL2 shows the existence of second copy of 
Tn3 immediately downstream of the defective Tn3 element, and so we 
have revisited the question of whether the opd cluster of pPDL2 
functions as an active transposon. Bioinformatic analysis identified 
Tn3-specific terminal repeats (Siguier et al. 2006) in pPDL2: two at 
the terminal ends of Tn3 and a third at the upstream region of trans- 
poson y4qE (Figure 5A). The genetic map of the DNA region between 
y4qE and Tn3 includes the previously reported opd gene clusters 
(Siddavattam et al. 2003), transposons Tn3 and y4qE, and the ligBA 
operon coding for the protocatechuate-4,5-dioxygenase alpha and beta 
subunits, along with a gene (orfl) coding for a LysR-type transcription 
regulator. The ligBA operon codes for protocatechuate-4,5-dioxygenase 
alpha and beta subunits showing 83% identity to the similar proteins 
of Xanthomonas campestris pv. campestris str. ATCC 33913 (NCBI 
reference no. NP_636196.1). While evaluating the functional status 
of lig operon, an expression plasmid (pSDP9) was constructed by 



placing it under the control of vector-specific transcriptional and 
translational signals. The resting cells prepared from induced cul- 
tures of E. coli BL21 (pDSP9) successfully generated ring-cleavage 
products from proctocatechuate, catechol, and benzenetriol indicat- 
ing the operon codes for a functional dioxygenase. As the cloning 
strategy translationally fuses ligA with vector sequences coding for 
a His-tag, the E.,coli BL21 (pDSP9) cells encode a LigA C -6His protein. 
Affinity purification yielded both LigA c _ 6His and LigB in approxi- 
mately equimolar concentration, revealing that these two function- 
ally related proteins form a complex (Figure 5A). Supporting the 
resting cell assay, the affinity-purified dioxygenase also generated 
ring-cleavage products from both protocatechuate and catechols 
(Figure 5D). In addition to LigB and LigA C -6His proteins, an addi- 
tional band with a size of 25 kD was enriched during affinity puri- 
fication (Figure 5C). It is not clear whether this 25 KDa protein is 
copurified due to specific interactions with Lig proteins. Protoca- 
techuate-4,5-dioxygenase is the key enzyme in the benzoate (hy- 
droxylation pathway) and 2,4-dichlorobenzoate degradation 
pathways (Adriaens et al. 1989). In addition to these two com- 
pounds, many aromatic compounds such as vanillate, isovanillate, 
phthalates, and benzoate derivatives are shown to be channelled 
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through the protocatechuate degradation pathway (Providenti et al. 
2006). When OP insecticides, like methyl parathion and parathion, are 
hydrolyzed by OPH, aromatic compounds such as p-nitroarophenols 
are generated (Sethunathan and Yoshida 1973). These compounds are 
found to be more toxic to the microflora than the parent compounds 
(US Environmental Protection Agency 2007). The existence of a ring- 
cleavage dioxygenase as part of the opd gene cluster allows channeling 
of these aromatic byproducts of OP insecticides into the TCA cycle, 
contributing to their complete mineralization. 

In addition to the lig operon, lysR and two ORFs showing homol- 
ogy to transporter proteins belonging to the major facilitator super- 
family were identified in the upstream and downstream region of the 
lig operon (Figure 5A). Transporters belonging to the major facilitator 
superfamily proteins are known to be involved in transport of aro- 
matic proteins (Pao et al. 1998). Hence, in summary, if the sequence 
found between the Tn3-specific terminal repeats is a transposable unit, 
it appears to have all necessary information for mineralization of OP 
insecticides like parathion and methyl parathion. 

The opd cluster of pPDL2 is an active transposon 

The transposition assay that we described in our previous study was 
repeated by replacing pMS3 with the pPDL2-K derivative pPDL2-KT in 
which opd is replaced with opdv.tet and hence the event of transposition 
can be monitored by the tetracycline-resistant phenotype. Plasmid 
pPDL2-KT requires pir protein for its replication and hence replicates 
only in E. coli pir-l 16 cells. When pPDL2-KT was electroporated into E. 
coli (pJQ210) cells, we found colonies resistant to sucrose and tetacy- 
cline, indicating disruption of the sacB gene of pJQ210 due to trans- 
position of opdv.tet. Further analysis of plasmids isolated from the 
sucrose/gentamycin/tetracyclin resistant colonies gave a typical restric- 
tion profile suggesting existence of opd cluster in the sacB region of 
pJQ210. The junctions of the opd cluster were then amplified using 
sacB- specific primers and primers designed for the region upstream 
of the predicted terminal repeats. The sequence of the PCR amplicons 
confirmed the existence of Tn3-specific terminal repeats, indicating that 
the opd cluster of pPDL2 is an active transposon (Figure 5B). 

DISCUSSION 

Horizontal gene transfer in bacteria is an important driving force that 
plays a critical role in generating genetic variation among bacteria and 
contributes to shaping of their genomes (Hacker and Carniel 2001; 
Jain et al. 1999, 2003; Preston et al. 1998; Skippington and Ragan 
2011). DNA acquisition by bacteria through various mobile genetic 
elements, such as plasmids, bacteriophages, IS elements, transposons, 
integrons, and conjugative transposons, is an established phenomenon 
(Dobrindt et al. 2004; Garcillan-Barcia et al. 2011; Wiedenbeck and 
Cohan 2011). Sequences associated with mobile genetic elements dif- 
fer from the host genome with respect to their GC composition and 
GC skew, serving as unique makers to assess lateral acquisition. Such 
delineation represents an island-like existence for mobile genetic ele- 
ments in the genome of the host (Juhas et al 2009). Depending on the 
nature of the genetic information present on these mobile genetic 
elements, they are designated as pathogenic islands (Dobrindt et al. 
2004; Schmidt and Hensel 2004), resistant islands (Larbig et al. 2002), 
or catabolic islands (Gaillard et al. 2006). Involvement of plasmids and 
catabolic transposons in lateral transfer of degradative genes is a recog- 
nized phenomenon (Topp et al. 2000; Tsuda et al. 1999; van der Meer 
et al. 2001). Recently genes required for 3-chlorobenzoate degradation 
have been shown to be part of an integrative conjugative element (ICE). 
This ICE element, found downstream of a glycyl tRNA gene, contains 
all the necessary genes for mineralization of 3-chlorobenzoate along 



with an int gene coding for a site-specific integrase and attachment 
sites (Gaillard et al. 2006). 

Plasmid pPDL2 is unique in number of ways. It contains a well- 
defined replication module with an oriV sequence and a rep A gene to 
support its replication. A mini-replicon generated by ligating the oriV, 
repA region of pPDL2 with an insertionally inactivated opd {opdv.tet) 
gene successfully replicated in Acinetobacter sp. DS002, suggesting its 
involvement in replication of pPDL2. The pPDL2 also has an integrase 
module, consisting of copG, int, and pgm genes (Figure 4A). The 
organization of the integrase module is similar to the one present in 
plasmid pUTl of S. japonicum UT26S, although pUTl has a single 
integrase module, whereas pPDL2 has two modules. The data pre- 
sented in the present study clearly show that pPDL2 can successfully 
integrate if a suitable attB site is available in the host chromosomal 
DNA. Plasmid pPDL2 contains an oriT sequence that is identical to 
the oriT sequence of pLBl of S. japonicum UT26. In support of this 
observation, the pPDL2 derivative pPDL2-K was mobilized into Aci- 
netobacter and E. coli cells in the presence of a helper strain E. coli 
HB101 (pRK2013). The data obtained in the present study clearly 
show formation of cointegrates between pPDL2-K and pSDP8 due 
to site-specific recombination. If a suitable attB site is present, such an 
event is possible between pPDL2-K and the E. coli chromosome. If 
this property of pPDL2 is taken together with its ability to mediate its 
own mobilization, designating pPDL2 as an IME appears to be logical. 
In general, lateral mobility of Sphingobium plasmids is only seen 
among Sphingomonadaceae members (Basta et al. 2004). However, 
the self-transmissible Sphingobium plasmid pLBl derivative pLBl-K 
has been successfully transmitted from S. japonicum to Sinorhizobium 
meliloti and Mesorhizobium loti (Miyazaki et al. 2006). Like pLBl, 
plasmid pPDL2 derivative pPDL2-K can be successfully mobilized 
from E. coli to Acinetobacter that is not a member of the Sphingomo- 
nadaceae. This property of pPDL2 facilitates lateral transfer of the opd 
cluster among a wide range of bacterial strains. When such lateral 
transfer takes place into a host where pPDL2 replication is not sup- 
ported, it can integrate into the genome of the recipient cells due to 
the existence of an integrase module. This additional feature appears 
to ensure retention of genetic information evolved for degradation of 
organophosphates, even in a host that does not permit replication of 
pPDL2. 

Furthermore, pPDL2 contains a catabolic transposon that contains 
information for hydrolyzing a structurally diverse group of organo- 
phosphates and then for degrading aromatic compounds generated 
due to OPH-mediated hydrolytic cleavage of OP insecticides. As 
shown in Figure 5A, the entire degradative module appears to be 
organized as a complex transposon having transposons y4qE and 
Tn3 at the flanking ends of degradative modules. Tn3-specific termi- 
nal repeats were found upstream and downstream of y4qE and Tn3, 
and consistent with this sequence information, we demonstrated 
a transposition event that inserted the entire degradative module into 
the sacB gene of a reporter plasmid. This event suggests the possibility 
of spreading opd as part of a large catabolic transposon. Surprisingly, 
most of the genetic information found on 37,317 kb pPDL2 appears to 
be related either to lateral gene transfer or degradation of organo- 
phosphates. It appears that the plasmid is tailor-made through a num- 
ber of recombination and transposition events for spreading opd 
information among soil bacteria. The presence of identical opd regions 
on plasmids pCMSl and pPDL2 and the existence of an active trans- 
posable element with opdA is in itself evidence for acquisition of opd 
through lateral gene transfer (Mulbry et al. 1987; Home et al. 2002, 
2003; Pandeeti et al. 2011). In the Sphingobiaceae family, only two 
members, S. fuliginis ATCC 27551 and Sphingomonas sp. JK1, have 
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opd genes (NCBI reference no. ACD85809). The existence organo- 
phosphates in the soil might have favored the opd-containing bacteria 
to acquire DNA that encodes transporters and aromatic degradation 
genes as it allows the complete utilization of organophosphates. 

Considering the experimental data presented in this study, we have 
searched genome sequences of Sphingobiaceae members to determine 
whether sequences that match pPDL2 exist in other genomes. Synteny 
maps were drawn to show the extent of similarity between pPDL2 and 
genome sequences of Sphingobiaceae members (S. japonicum and 
S. wittichii). Plasmid pUTl of S. japonicum has organizational simi- 
larity with pPDL2 in the region containing replication and integrase 
genes (see Figure S2). If synteny maps are taken into consideration 
along with GC content and GC skew, it clearly suggests that plasmid 
pPDL2 evolved by a number of transposition and rare recombination 
events to integrate DNA from different sources. Although pPDL2 has 
all the necessary features for spreading opd information, due to the 
lack of tra genes, it depends on helper plasmids for achieving lateral 
transfer. However, as Sphingobium fuliginis ATCC 27551 has three 
more plasmids, complementation of such a function by these plasmids 
cannot be ruled out. 
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